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ABSTRACT 

Intermediate mass stellar evolution tracks from the main sequence to the tip 
of the AGB for five initial masses (2 to 6 M©) and metallicity Z=0.0001 have been 
computed. The detailed ID structure and evolution models include exponential 
overshooting, mass loss and a detailed nucleosynthesis network with updated 
nuclear reaction rates. The network includes a two-particle heavy neutron sink 
for approximating neutron density in the He-shell flash. It is shown how the 
neutron-capture nucleosynthesis is important in models of very low metallicity for 
the formation of light neutron-heavy species, like sodium or the heavy neon and 
magnesium isotopes. The models have high resolution, as required for modeling 
the third dredge-up. All sequences have been followed from the pre-main sequence 
to the end of the AGB when all envelope mass is lost. Detailed structural and 
chemical model properties as well as yields are presented. This set of stellar 
models is based on standard assumptions and updated input physics. It can 
be confronted with observations of extremely-metal poor stars and may be used 
to assess the role of AGB stars in the origin of abundance anomalies of some 
Globular Cluster members of correspondingly low metallicity. 



Subject headings: stars: AGB — stars: evolution — nuclear reactions, nucle- 
osynthesis, abundances — stars: interiors 
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1. Introduction 

The chemical production of low- and intermediate mass stars during their asymptotic 
giant branch phase of evolution is an important contribution to the galactic chemical evo- 
lution. For example, half of all elements heavier than mass number A = 90 are made in 
low mass AGB stars by the s-process. AGB stars are responsible for a significant fraction 
of the carbon in the solar system abundance mix. AGB stars are the favored candidate 
for the main primary nitrogen source at lower metaUicity (Pettini et al. 2002). The dis- 
crepancy of observations and galactic chemical evolution predictions based on massive star 
yields is in particular evident for some neutron heavy isotopes of certain elements (like Mg 
or Na, Timmes et al. 1995). This indicates the role of very low metallicity AGB stars to the 
early galactic chemical evolution. Equally important are the informations obtained from an 
increasing number of spectroscopic observations of extremely metal-poor (EMP) stars. 

In the past yield predictions for a range of initial masses and metallicities have been de- 
rived from so-called synthetic AGB models (Renzini & Voli 1981; Marigo et al. 1996; van den 
Hoek & Groenewegen 1997). In these models the evolution is described by constructing fit- 
ting formulas to certain quantities of full stellar evolution models that solve the full set of 
stellar structure and nuclear energy generation. The most recent improvements include a 
full envelope integration for the hot bottom burning phase of massive AGB stars (Marigo 
1998). The synthetic AGB model approach is very well justified under certain circumstances 
and for specific purposes. The computation of AGB stellar models requires high numerical 
resolution and significant computing time. In the past large grids in mass and metallicity 
with several ten- to onehundredthousand models for each sequence have not been feasible 
due to their computationally demanding nature. This is changing now and as an example, 
Karakas et al. (2002) have presented a comprehensive grid of AGB calculations to study 
differentially the dependence of DUP on mass and metaUicity. Still, full stellar evolution 
models in the past were not able to reproduce, for example, the third dredge-up (DUP) in a 
way that agrees with observations. New results from nuclear physics may help to resolve this 
issue in the future (Herwig & Austin 2004). Synthetic models on the other hand can param- 
eterize the third DUP in some way, and calibrate its efficiency with observed properties of 
AGB stars, for example using the C-star luminosity function of the Magellanic Clouds. The 
synthetic model predictions arc therefore internally consistent with observations. In addition 
these models arc useful because they condense the information contained in the C-star LF, 
or in C-rich star and 0-rich star number counts in the Magellanic Clouds into more simple 
constraints, which detailed structure and evolution models have to reproduce. For example, 
the synthetic models by Marigo et al. (1996) require that efficient dredge- up must take place 
at core masses as low as 0.58 Mq, a value which has even been revised downward by a more 
recent update (Marigo et al. 1999). 
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In order to model yields and the abundance evolution for EMP stars full stellar models 
appear to be in particular important. We are just beginning to identify the peculiarities 
of nucleosynthesis in EMP AGB stars (or at Z=0 should they have existed, Chiefii et al. 
2001). Synthetic models rely both on accurate structure and evolution models to derive the 
fitting formulas, as well as on a good calibration. Both are presently not well established 
at metallicities of Z < 0.001. Models of massive AGB stars at very low metallicity for 
example, show behaviors not seen in models of larger metallicity, like H-burning during the 
DUP or DUP reaching below the He-shell (Herwig 2004). It is therefore not appropriate to 
simply extrapolate synthetic models that have been calibrated at moderate metal deficiency 
to very low metallicities. In order to improve the situation more complete stellar evolution 
calculations are needed. Here, the structure and chemical evolution for a set of five tracks 
with masses from 2 Mq to 6 Mq at metallicity Z — 0.0001 are presented. 

The application of these models may include comparison to observations of EMP stars 
or tests of the possible role IMS (intermediate mass stars) in the chemical evolution of low- 
metallicity Globular Clusters. Although the grid is far from complete in mass and metallicity 
the calculations may be useful for Galactic Chemical Evolution models. Since the structural 
evolution is given as well, the models may be used to investigate the s-process. 

In Paper I (Herwig 2004) the evolution of massive AGB stars at very low metallicity has 
been described in detail, in particular the interplay of hot bottom burning and the hot DUP. 
In another previous paper we have specifically addressed the possible implications of the 
more massive cases at this metallicity for the star-to-star abundance variations observed in 
globular clusters (Denissenkov & Herwig 2003). Here the detailed thermodynamic, structural 
and abundance evolution of a homogeneous set of intermediate mass EMP stellar models is 
presented. In § 2 the code changes compared to the version used in Paper I are described. 
The model results are given in § 3. The conclusions are presented in § 4. 

2. Physical input, stellar evolution code and model calculations 

The stellar evolution code and all model assumptions are the same as in Paper 1, with 
the following exceptions. 

2.1. Mass loss 

Mass loss for AGB stars of very low metallicity is neither observationally nor theoret- 
ically well constrained, and one has to rely on extrapolation from larger metallicities and 
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rather qualitative and indirect considerations. For a mass loss formula both the functional 
dependence of mass loss on stellar parameters as well as the absolute calibration are needed. 

Mass loss rates of Mira-like pulsating stars can be calculated by modehng the dynamical 
structure and evolution of their atmospheres (Bowen 1988). For solar-hke metallicity Blocker 
(1995) has extracted the luminosity dependence of the Bowen (1988) models and finds 

logM(x3.71ogL. (1) 

The physical input like the opacities as well as the properties of dust formation of the 
dynamical atmosphere models depend on the metallicities. However, no such models are 
available for low or very low Z. Here I assume not only that Eq. (1) can be applied at 
Z = 10""^, but also that it approximately applies to changes in L due to changes in Z too. 
This assumption may in fact not be entirely foolish. The DUP is more efficient in EMP 
AGB stars (Herwig 2004, and § 3) and the CNO elements collectively reach close to the solar 
value abundances in the envelope. A comparison of the Rosseland mean opacity for a 3 Mq 
TP- AGB model at Z = 0.02 and Z = 10"^ yields a difference of typically less than 10% for 
logT < 4.6. 

If the dependence Eq. (1) is applied to luminosity variation due to varying metallicity, 
then the variation of stellar parameters with metallicity introduces a significant variation of 
the mass loss. In fact due to this effect one obtains a larger mass loss at lower Z because 
the luminosity is larger at lower Z. The average luminosity and temperature of a 2 Mq and a 
5 Mq for a range of metallicities is shown in Fig. 1. From these stellar structure and evolution 
results one determines approximately 

logL oc -0.21ogZ . (2) 

Together with Eq. (1) this implies that M{Z = 10"^) ~ % ■ M{Z = 0.02) with rji ~ 50. 

van Loon (2000) has attempted to derive a mass loss-metallicity relation based on ob- 
servations of Magellanic Cloud giants, and suggests that 

logM oc0.31ogZ. (3) 

This exponent has a very large uncertainty, mainly due to small-number statistics. Due to its 
observational nature Eq. (3) contains both the possible dependence of Eq. (1) on metaUicity 
and the dependence of stellar parameters on metallicity (Eq. 2). Tentatively extrapolating 
the relation of van Loon (2000) to Z = 10~^ the absolute mass loss rates ai Z — 0.02 and 
Z = 10"^ are related by M{Z = 10"^) ~ r]2 ■ M{Z = 0.02) with rj2 ~ 0.2. 

One can then obtain a cahbration for the Bowen-Bloecker mass loss rate (Eq. 1) at 
Z — 10"^ by 77b = V2/V1 which gives 77B = 0.004 for the exponent 0.3 of Eq. (3). However, it 
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has been mentioned that this exponent is subject to large uncertainty. In particular, both 
result for M stars and for C stars are compatible with logM oc — 0.3 log Z which implies 
rj2 ~ 5. Therefore the largest mass loss at Z = 10"'^ within this framework is t/b = 0.1 and 
this is the adopted value for the present study. In this way these calculations give a lower 
limit to the envelope enrichment and the corresponding yields. A lower mass loss rate would 
lead to a longer TP-AGB duration with more TPs and subsequent DUP events. 

2.2. Overshooting 

The treatment of convective boundaries, both at the bottom of the envelope convection 
as well as at the bottom of the pulse-driven convection zone (PDCZ) during the Hc-shcll flash 
effects stellar models in many ways, including the dredge-up predictions. The consequences 
and implications of various treatments of mixing at these convective boundaries has been 
explored in detail (Herwig et al. 1997; Mowlavi 1999; Herwig 2000; Langer et al. 1999; Lugaro 
et al. 2003; Herwig et al. 2003; Herwig 2004). Models which include an adjustable amount 
of depth- and time-dependent amount of overshooting at all convective boundaries show 
efficient dredge-up at low core-masses. 

The same concept of exponential overshooting as in Paper I is applied. For the TP- 
AGB models the overshoot efficiency at the bottom of the PDCZ is set to /pdcz = 0.008 
while the efficiency at the bottom of the convective envelope is set to /ce = 0.016 at all 
times. Models with larger /pdcz have higher temperature for the ^^Ne neutron source. 
Lugaro et al. (2003) found that if /pdcz is too large some s-process branchings produce 
isotopic ratios that are irreconcilable with laboratory measurements in pre-solar SiC grains. 
For example, the branching at ^^Zr is sensitive to the temperature at the bottom of the 
He-shell flash convection zone. If /pdcz is too large the temperature is too large and the 
predicted ^^Zr/^^Zr ratio is larger than measured in the grains. Lugaro et al. (2003) have 
concluded that /pdcz = 0.016 may be too large. Unfortunately the neutron cross section of 
the unstable isotope ^^Zr is very uncertain. The current range of estimates differ within a 
factor of four (Kawano, priv. com., Bao et al. 2000a; Jorissen & Goriely 2001), and a more 
conclusive analysis has to await a better experimental determination of this cross section. 
Other nuclear reaction rates as studied by Herwig & Austin (2004) may also be important. 
In any case the concept of overshooting at the bottom of the PDCZ plays an important role 
in currently favored models of H-dcficicnt central stars of planetary nebulae (of spectral type 
PC 1159 and [WC]-CSPN, Koesterke & Hamann 1997; Dreizler et al. 1996; Werner & Herwig 
2005). 

The rather small value for /ce = 0.016 can not reproduce the partial mixing needed 
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for the ■'^^C neutron source in low mass TP-AGB stars (Goriely & Mowlavi 2000; Herwig 
et al. 2003). One can estimate the mass of the s- process layer resulting from partial mixing 
that reproduces the observed stellar s-process overabundances. The stellar model must 
accomplish such an enrichment within a number of dredge-up events that can be constrained 
by observations. This number can be derived from model comparison with the observed 
C-star luminosity function and the observed ratio of C-stars to 0-rich AGB stars (Marigo 
et al. 1999). In particular newer synthetic models including the effect of C/0 ratio dependent 
molecular opacities (Marigo 2002) indicate that typical low-mass AGB stars that produce 
the s-process elements may in fact experience only five to at most ten thermal pulses and 
subsequent dredge-up events during the C-rich phase. 

In the context of evaluating the properties of rotating AGB star models Herwig et al. 
(2003) estimated that the partial mixing zone that hosts the s-process nucleosynthesis should 
have a mass Mp > 7 • lO^^M©, in agreement with more detailed calculations. Exponential 
overshooting as applied here with an efficiency /ce = 0.016 leads to a partial mixing zone 
with a mass ^ 1O~^M0 only (Herwig 2000; Herwig et al. 2003). In a detailed analysis 
Lugaro et al. (2003) have assumed that /ce = 0.128 and obtained overall overabundance 
factors close to what is observed. There are indications that /ce may be even larger in the 
framework of the s-proccss partial-mixing concept, at least for solar-metallicity cases. 

However, in Paper I it was shown that in models of extremely low metallicity even a very 
small amount of exponential overshooting leads to hot DUP with potentially very significant 
implications to the overall evolution. In these stars the ^^C rich intcrshcll material beneath 
the proton-rich matter is hotter than in cases with higher metallicity. If protons are forced 
into the ^^C rich core during the dredge-up phase of the TP cycle they are burnt vigorously. 
The additional H-burning luminosity can further drive convective instability and cause a 
corrosive, convective H-shell that thus penetrates much deeper into the core than without 
this additional H-burning. This hot DUP and the resulting corrosive burning is in fact 
conceptually related to the conductive propagation of nuclear fiames considered in ONeMg 
cores in the context of SN Type la models (Timmes et al. 1994). 

In Paper I calculations of a 5 dredge-up phase with /ce = 0.03 - about twice as 
large as the value assumed here - prompted ongoing dredge-up that would have terminated 
the AGB in approximately 2000 yr (based on some uncertain assumptions on mass loss for 
such configurations). The situation for lower- mass cases may be less dramatic, but still 
hot DUP may alter the conditions for the formation of a partial mixing zone at extremely 
low metallicity, as enhanced H-burning during this phase would result from any efficient 
convection-driven extra-mixing process. 

While it is possible that exponential extra-mixing as assumed in the diffusive over- 
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shooting approach is responsible for the s-process partial mixing zone at solar metallicity, 
the situation is rather unclear at extremely low metallicity. The situation can be improved 
by investigating how the structure (i.e. mainly the DUP depth) depends on the interplay of 
extra-mixing and hot DUP, and how extra-mixing itself is effected by nuclear burning in the 
overshooting layer. In these calculations a rather small efficiency of overshooting is assumed. 
This way most of the complications caused by hot DUP are avoided. However, the nuclear 
production predictions do not include the uncertain effect of a more extended ^^C pocket. 



2.3. Nuclear network 

The updated nuclear reaction network includes all relevant charged particle reactions, 
/5-dccays and neutron-capture reactions. If available the NACRE reaction rates arc used 
(Table 1). The nuclear network solution is based on rates from linear interpolations in the 
250 grid point logarithmic tables generated by the WWW database tool NETGEN ( Jorissen 
& Goriely 2001). For the calculation of the energy generation in the stellar structure solution 
a smaller network containing the dominating CNO cycle, pp-chain, and He-burning reactions 
has been used. The reaction rates have been calculated from the fitting formulas to the 
NACRE reactions (Angulo et al. 1999), which provide smooth T-derivatives. Potentially this 
dual approach may introduce an inconsistency because a small amount of energy released 
by trace elements may be ignored, and because the fitting formulas may disagree with the 
tabulated values by a few percent. However, it turns out that the fitting formula (despite 
the published fitting errors) represent the tabulated values very well for the most important 
H- and He-burning reactions in the relevant temperature range. 

All important light n-capture reactions as well as the several Fe-group species capturing 
neutrons are included. A light and a heavy n-sink take care of the remaining n-capture species 
as in Herwig et al. (2003). The two-particle sink approximating the trans- iron species allows 
a rough estimate on the number of neutrons captured per iron seed particle (ricap) which is 
an important diagnostic tool for the s-process. For the sink-cross sections the same value as 
in Herwig et al. (2003) is taken. This cross section depends on the heavy elements abundance 
distribution, and therefore a cross section derived for (near-)solar metallicity is not correct 
at very low metallicity. However, a detailed study of the s-process at very low metallicity 
which is required to estimate an appropriate sink-cross section is not available at this time. 
In addition the exact value of the heavy sink-cross section is less important at very low 
metallicity because the role of heavy species for the determination of the neutron density is 
smaller when large amounts of primary light n-capture species, like ^^Ne, are present. 
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3. 



Results 



3.1. 



Structural evolution 



The average structural properties of the model grid are summarized in Table 2. It shows 
how the stellar parameters as well as evolutionary times change systematically as a function 
of initial mass. 

The surface temperature of the models as a function of TP-phase and mass changes 
only little, and all TP-AGB models are in the range Teff = 4000 — 4350 K. All sequences 
show very efficient third dredge-up. Therefore the core masses are on average almost con- 
stant during the TP-AGB evolution. These full stellar evolution models do not support the 
scenario of low- mass supcrnovae in the early Galactic halo proposed by Zijlstra (2004). An 
important assumptions in this proposal is that the cores of IMS stars are able to grow to 
the Chandrasekhar mass because the mass loss is very low at low mctallicity. The models 
presented here can not rule out this proposal, but it is fair to say that they do at least 
raise some serious doubts. Deep dredge-up events prevent the core from growing. They also 
cause substantial pollution of the surface layers with C, N and O, all of which are important 
sources of opacity in the outer stellar layers. Thus the assumption that mass loss is in fact 
neghgible in these stars may not be correct. Finally, the hot DUP (previous section and 
Paper I) is more effective at higher core mass. Thus, from these calculations I have more 
reasons to doubt rather than believe this idea. 

The stellar luminosity and radius evolve only httle with TP number and during a TP 
cycle. Averaged values are given in Table 2. These TP-AGB models with initial masses of 
2 ... 6 M0 have core masses in the range 0.62 ... 1.05 Mq. The luminosities follow approxi- 
mately an exponential core-mass luminosity relation: 



This is shown in Fig. 2, where in addition two data points for TP-AGB stars (last TP cycle) 
of zero- metalli city by Chieffi et al. (2001) arc shown. Our models show sharp luminosity and 
radius peaks during the first deep dredge-up events (see Paper I for details). These peaks 
may be important if the binary evolution and of EMP stars (in the GH-star scenario) is 
considered. It may be that interaction or enhanced tidal synchronization occurs preferentially 
during these spikes close to the time of the He-shell flash. 

The core mass is defined here as the mass within which the hydrogen mass fraction is 
Xn < 0.37. The core mass grows because of H-shell burning. The recurrent third dredge-up 
events decrease the core mass. The result is a characteristic saw-tooth shaped curve of the 
core-mass as a function of time. In Fig. 3 the core mass evolution is shown for all sequences. 



log L/ Lq = 1.7414 X MJ Mq + 2.8799 



(4) 
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The DUP parameter is defined as A = AMh/AMdup, where AMh is the core mass growth 
due to H-burning during the interpulse phase. In all sequences A ~ 1 for the majority of 
TPs. Models with larger core mass have a shorter interpulse phase and smaller intershell 
layers (see Table 2), as it is the case at larger metallicity. For the mass loss law chosen here 
the more massive models experience more DUP events during which less processed material 
is mixed into the surface. This in addition to the larger dilution factor of the massive AGB 
stars due to their larger envelope mass, leads to the prediction that more massive AGB stars 
will generally have smaller overproduction factors, for example for the overall enrichment in 
the sum of CNO elements. 

In addition to the mixing properties the temperatures at the nuclear production sites 
determine the chemical enrichment of the stars. In Table 3 the temperatures in the H-shell, 
the He-shell as well as the temperature at the bottom of the envelope convection and the 
bottom of the PDCZ are given. The temperature in the He-shell and the temperature at 
the bottom of the envelope convection zone vary substantially and the values given in the 
table should be considered typical numbers taken approximately when half of the interpulse 
phase has passed. The last line for each sequence gives the core mass and envelope mass 
when the computations have been stopped. The subsequent evolution would lead into the 
post- AGB and central star of planetary nebula phase of evolution. The sequences have been 
stopped before this numerically difficult transition phase away from the AGB starts, but 
after enough envelope mass has been lost for reliable yield predictions. 

The density at the bottom of the PDCZ ranges from 4.0 x lO^g/cm"^ to 4.8 x lO^g/cm^ 
for the 2 Mq sequence (E82) while slightly lower densities are encountered at higher core 
masses (2 x 10^ - 1.8 x lO^g/cm^ for the 6Mq sequence E86). 

3.2. Chemical evolution and yields 

In the context of Galactic Chemical Evolution different species evolve differently as a 
function of metallicity (time) and depending on their particular location. In order to choose 
initial abundances for the stellar evolution models one can either adopt theoretical initial 
abundances from pubhshed galactic chemical evolution models (Timmes et al. 1995), or the 
initial composition is tailored to reproduce observed abundances at [Fe/H]= —2.3. Unfortu- 
nately, both approaches lead to different results for some important species (hke N or Na), 
while the theoretical models can not be checked in other cases (like ^^Ne). Fortunately the 
yields of many important species of TP-AGB stars are dominated by primary nucleosynthe- 
sis production. Moderate uncertainties in the initial abundances will have only little effect 
on the results. This is certainly true for the CNO elements. Therefore, the initial composi- 
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tion is set to metallicity scaled solar abundance distribution. The largest difference between 
the initial metallicity at Z = 10~^ and a solar-scaled initial composition would concern the 
light elements deuterium and ^He. These should be overabundant at very low metallicity 
compared to the solar-scaled values due to big bang nucleosynthesis (Walker et al. 1991). In 
this study the light elements lithium and ^He are not considered. 

IMS experience several kinds of dredge-up events that generate overabundances in their 
envelopes. At extremely low metallicity the second dredge-up after the end of core He- 
burning plays a more important role than the first dredge-up after the core H-burning. In 
fact in this model set only the 2 Mq sequence has a first dredge-up that is deeper than the 
second dredge-up. The most massive cases do not even have a first dredge-up and Red Giant 
Branch evolution phase (Girardi et al. 1996). He-core burning sets in before the model star 
reaches the giant branch and temporarily reverses that evolution (Fig. 4). After that they 
evolve directly into early AGB stars and the envelope becomes convectively unstable for the 
first time. Although the repeated third dredge-up events are mainly responsible for the total 
chemical envelope enrichment with nucleosynthesis products it is nevertheless important to 
document the pollution during the evolution prior to the thermal pulse AGB. 

The overabundances just before the first dredge-up event of the thermal pulse AGB 
evolution (Fig. 5) reflect mostly abundance changes that are secondary in nature. The ■''^C 
abundance is in all but the 6 M0 case depleted by mixing the partially CN cycled envelope. 
Accordingly the -"^"^N abundance is enhanced. In the 6 Mq sequence the Hc-shcU is hot and 
broad enough after the end of He-core burning that the descending envelope convection 
(second dredge-up) engulfs a small amount of primary ^^C. This explains also other small 
CNO abundance differences between the 6 Mq and the less massive cases. The Ne-Al isotopes 
are effected by secondary p-capture nucleosynthesis as well. This is more so the case for larger 
masses. In particular ^^Na is produced from p-captures on Ne isotopes. For this element 
the production prior to the TP- AGB is significant compared to the later production during 
the TP phase if a solar-scaled ^^Ne initial abundance is assumed. However, this may not 
be correct since standard models of massive stars (the yields of which may dominate the 
initial abundance distribution of EMP low- and intermediate mass stars) predict a lower 
than solar-scaled ^^Ne abundance. 

To compare the overall abundance evolution as a function of mass the envelope mass 
averaged mass fractions have been computed for each sequence and each species i using 



where Mi and Mf are the stellar mass at the beginning and the end of the AGB evolution 
phase, and X is the surface mass fraction at a given time corresponding to the stellar mass. 




(5) 
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These averaged abundances in the matter returned to the interstellar medium together with 
the adopted initial abundance distribution is given in Tabic 4. These averaged overabun- 
dances should be useful for comparison with trends in observed abundances of metal poor 
globular clusters or extremely metal-poor halo stars. Note that the squared bracket notation 
[X/Fe] is very well approximated by logXav/Xini because the initial abundance is solar- 
scaled. The average overabundances are displayed for all cases in Fig. 6. Yields according 



are given in Table 5^. 

The evolution of surface abundances and stellar parameters for all tracks is provided 
in abbreviated online tables (see Table 6 for a sample) . A graphical representation is given 
in Fig. 7 to 10. For the following discussion two abbreviations are useful: LC refers to the 
low-mass cases 2 and 3 which do not experience hot-bottom burning, while HC refers to 
the 4-6 Mq cases in which HBB is efficient. 

Helium ^He is brought into the envelope by all dredge-up episodes. In addition HBB in 
HCs transforms H into helium. HCs are more efficient ^He producers than LCs, both in terms 
of the yields as well as the average abundance in the ejecta. The relative importance of the 
different production mechanisms is a function of initial mass. HCs generate the vast majority 
of "^He during the second dredge-up. Although HBB is most efficient in the most massive 
cases the 6 sequence (E0086), for example, generates only a few percent of the total "^He 
overabundance during the TP-AGB phase. Contrary, the 2 case (E0082) produces about 
two thirds of its ^He overabundance during the TP-AGB by third dredge-up. The mass 
dependence of the ^He yields suggests that for a Salpeter-like intial mass function both HCs 
and LCs are equally important producers of ^He. 

Carbon ^^C is produced in the He-shell by the triple-o; reaction and dredged-up after the 
TP. The LCs show larger ^^C overabundances and also larger yields then the HCs. The 
deeper dredge-up (in mass, not in dredge-up parameter A) and the smaller dilution factor 
due to a smaller envelope mass outweigh the fact that LCs have less TPs than HCs. ^^C is 
not produced or destroyed in the 2 case which does not have HBB. The matter which is 
dredged-up from the intershell is void of ^^C beacuse of the large ^'^C(q;, n)^^0 cross-section. 
For this mass the overabundance is entirely secondary and stems from the pre-TP-AGB 



Table 5 available online. 



to 




(6) 
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evolution. For larger masses HBB becomes more and more important and so does the ^^C 
overabundance, which is primary in these cases and exceeds the enrichment level obtained 
by the second dredge-up. For M > 4 Mq the HBB is efficient so that the ^^C overabundance 
exceeds the ^^C overabundance by about a factor of 10. Both low-mass stars as well as 
massive stars produce carbon with much larger isotopic ratio than this (Woosley et al. 
2002). A nuclear production site like these intermediate mass stars with a small ^^C/^"^C 
ratio is therefore needed to account for example for the solar ratio of ^^C/^^C ~ 89. 

Nitrogen The LCs maintain their isotopic nitrogen enrichment from the ffist and second 
DUP. The nitrogen overabundance in these heavily C-enriched model stars is therefore of 
secondary origin. The upper limit of the nitrogen abundance which can be expected in stars 
polluted by the LCs is given by the initial metaUicity. 

Large primary production of ^^N are resulting from combined efficient DUP and HBB 
for the HCs only. The 5 and 6Mq models even show a small ^^N enrichment. The ^"'N/^^N 
ratio is determined by the CN cycle equilibrium value. If DUP and HBB are efficient then 
the ^^N overproduction is so large that some primary ^^N is produced as well. However, this 
does not address the problem of the nuclear production site of the solar ^^N abundance. In 
the absence of detailed models earlier studies have speculated that ^^N may be produced 
in IMS (Timmes et al. 1995). This possibility can be ruled out at least for models which 
employ rather standard physical assumptions. 

Finally note that the carbon to nitrogen ratio in the ejecta of IMS stars is almost 
bimodal. The LCs have large [C/N] while the HCs have small ratios. With standard as- 
sumptions on mixing the mass transition between high and low [C/N] is very sharp. It is 
therefore not clear that the simultaneous C and N enrichment observed in many EMP CH- 
stars can be explained by mass transfer from companion stars that have evolved like any of 
the models presented here. 

Oxygen ^^O is enhanced in the envelopes of EMP AGB stars by third dredge-up. The 
pattern is the same as in the case ^^C: the lower mass cases bring cumulatively more intershell 
material into the envelope compared to the higher mass cases. The origin of ^^O is the 
^^C(q;,7)^^0 reaction transforming primary ^^C from the triple-a reaction. ^^O is mainly 
produced in the lower part of the He-shell where the ^He abundance is declining. During 
the He-shell flash the triple- a reaction dominates. The amount of ^^O available in the 
convectively mixed intershell depends on the treatment of convective boundaries. With our 
assumption of exponential overshoot at the bottom of the PDCZ an ^^O abundance in this 
layer of about 6% by mass is found, depending somewhat on TP number and mass. As shown 
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in Paper I the prediction of significant dredge-up of ^^O in EMP AGB stars is not a result of 

assuming overshooting at the bottom of the PDCZ. Even without this kind of extra-mixing 
the ^^O abundance in the PDCZ in EMP stars is of the order 1 — 2% by mass. This primary 
oxygen in the intershell is in any case several orders of magnitude more abundant than the 
initial ^^O abundance in the envelope. The importance of ^^O dredge-up in AGB stars of the 
lowest AGB stars is also found in the calculation for metallicity Z=0 by Siess et al. (2002). 
For example at the end of the evolution of their 1.5 M0 model which does not experience 
HBB the surface ^^C-mass fraction is 6.9 ■ 10"=^ and that of ^^O is 2.4 ■ 10"=^ 

The ^^O abundance is depleted rather rapidly for the HCs during HBB and not re- 
plenished because the ^^0(p, a^^N reaction exceeds the ^''0(p, aY^F reaction by almost two 
orders of magnitude. 

Neon A small primary production of ^°Ne can be noted for all masses (Fig. 9). This is just 
a curiosity as the nuclear origin of ^°Ne is carbon burning in supernova Type 11. However, 
it is noted that the small increase of ^'^Ne observed in these models is due to two reactions: 
1^0(0;, 7)2°Ne and to a lesser degree ^^0(n, 7)^^0(0;, n)20Ne. 

In the LCs the primary ^^Ne can increase the elemental neon abundance significantly. 
Both heavy neon isotopes are destroyed by HBB in the HC cases. During the final pulse cycles 
HBB becomes less efficient. Only then is the DUP of ^^Ne and ^^Ne more efficient than the 
interpulse destruction. In the He-shell fiash ^-"^Ne is mainly made by ^®0(n, 7)^'''0(q;, 7)^-^Ne 
and to a small extent by ^°Ne(n, 7)^^Ne. In the LCs in particular the ^^Ne yields are signifi- 
cant. It is important to note that even if the ^°Ne(n, 7)^^Ne dominates ^^Ne is still primary, 
because ^°Ne is produced in a primary mode, even though this is not important for the ^^Ne 
yield at this metallicity. 

Sodium Sodium is produced by both LCs and HCs by different processes. In HCs ^^Ne is 
dredged up from the intershell and then transformed into ^^Na during HBB. If, like in the 
6 Mq case, HBB is very efficient then ^"^Na is first produced and destroyed in the second half 
of the interpulse phase. However in all other cases significant amounts of ^^Na are produced. 
In the LCs the main source of ^^Na is ^^Ne(n, 7)^^Na. 

Magnesium The different nuclear production sites during H- and He-shell burning in TP- 
AGB stars is discussed in all detail in Karakas & Lattanzio (2003). Although they have 
computed models for somewhat higher metallicity all their basic mechanisms apply equally 
here. ^^Mg is only significantly changed if HBB is very efficient, like in the E0086, GMq 
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case (Fig. 10). Then both ^^Mg as well as ^^Mg are produced. This production mode of the 
neutron heavy isotopes of magnesium has a distinct signature: the ^^Mg abundance is larger 
than the ^^Mg abundance (Karakas & Lattanzio 2003; Denissenkov & Herwig 2003). Even 
taking into account nuclear reaction rate uncertainties does not change this finding. 

There is a second production site for the neutron heavy magnesium isotopes. In the 
He-shell flash convection zone the two a-captures on ^^Ne arc almost equally important: 
^^Ne(«, n)^^Mg and ^^Ne(Q;, 7)^^Mg. However, the temperature dependence of these two 
reactions is different. For Tg > 2.8 the (a, n) reaction exceeds the («,7) reaction, while 
for lower temperatures ^^Ne is more effectively transformed into ^^Mg. This dividing tem- 
perature coincides roughly with the activation temperature of the ^^Ne n-source reaction. 
Therefore, for low temperatures in the He-shell flash convection zone (as can be found in low 
mass models) one expects only little neutron-heavy magnesium production favoring ^^Mg. 
However, at higher He-shell flash temperatures should ^^Mg dominate the heavy Mg isotope 
production. 

There is a component missing in this picture, which appears to be important in EMP 
stars. ^^Mg has a very high neutron cross section, about 600 times larger than that of ^^C, 
about 200 times larger than ^^O, and still about twice as large as the ^^Al n-capture cross 
section. The relative number of neutrons actually available for ^^Mg depends not only on 
the relative value of the cross sections but equally on the abundance of competing neutron 
capturing species. Assuming, for example, that 5% by mass of a 2% mass abundance of ^^Ne 
is transformed into ^^Mg implies a mass abundance for this isotope of 10-"^. Thus, for this 
example (< av >(n,7) ^)25Mg exceeds (< av >(n,7) ^)i2c- Thus, neutrons are captured by 
^^Mg in significant numbers if ^^Ne releases neutrons. In order to check this picture I have 
conducted some one-zone nuclear network calculations for typical conditions found in the 
He-shell flash convection zone. These calculations confirm that neutron captures on ^^Mg 
are important for the evolution of the heavy magnesium isotopes in EMP AGB stars. For 
temperatures where ^^Ne releases neutrons (in the test calculations Tg = 3.2) X(^^Mg) > 
X{^^Mg) if all neutron capture reactions are considered, but X(^^Mg) < X(^^Mg) if the 
^^Mg(n, 7)^^Mg reaction is switched off, as expected from the ratio of two ^^Ne a captures. 

It should also be pointed out that this ^'''Mg and ^^Mg from the Hc-shell is of pri- 
mary origin while ^^Mg and ^^Mg from HBB is mainly secondary. There it is produced at 
the expense of the initial ^^Mg abundance, that has its primary origin during carbon and 
subsequent neon burning in Type II supernova. 

These models show that in fact EMP AGB stars can - within the uncertainties of mass 
loss in particular - produce primary neutron heavy magnesium isotopes with a ^^Mg/^^Mg 
ratio exceeding unity. This is in particular the case when mass loss is efficient and prevents 
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HBB from dominating the nucleosynthesis of the neutron-heavy Mg isotopes. This result 
should be relevant for the interpretation of recent isotopic magnesium abundance determi- 
nations ( Yong et al. 2003) . 

More comparison with other calculations Ventura et al. (2002) have presented IMS 
yields from complete stellar evolution calculations down to a metallicity of Z=0.0002. They 
focus on the evolution of helium and the CNO elements and have not included any heavier 
elements. Some of their input physics assumptions are different than those adopted in this 
study. These include in particular their choices on the treatment of convection. They use 
the Full-Spectrum Turbulence theory (Canuto & Mazzitelli 1991) which tends to cause more 
efficient HBB compared to the standard MLT. More importantly, they chose to apply no 
overshooting or any other means of extra-mixing during the TP-AGB phase. Even at very 
low metallicities and the corresponding high core masses their third dredge-up is as small as 
A = 0.2 . . . 0.4. As a result their average masses in the ejected material of species hke ^He, 
^^C or ^^O are significantly smaller than in this study (Fig. 11). A probably small part of 
this difference may be attributed to the metallicity difference of a factor of two between the 
two sets. However, while the two sets of models are quantitatively different the abundance 
trends with mass are the same. 



4. Conclusion 

This paper describes the properties of a grid of intermediate mass stellar evolution 
models with low metallicity (Z=0.0001). The computations include an extensive nuclear 
reaction network which allow predictions of yields and surface abundance evolution for a 
large number of species. The model also includes neutron capture reactions which appear 
to be important for the correct modeling of light neutron heavy isotopes. A rather small 
amount of exponential overshooting has been included. No other non-standard effects (hke 
rotation or magnetic fields) have been considered. Thus, these are rather standard models 
with updated input physics. 

Confronting models and observations at very low metaUicity is especially interesting as 
Galactic Chemical Evolution had less time to cover the tracks of individual nuclear pro- 
duction events. In this lies the importance of studying stars of the lowest metallicity, both 
theoretically and observationally. Confronting these standard models with observations of 
EMP stars will help to identify what additional processes are important. 

The amount of exponential overshoot in these models is too small to generate a ^^C 
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pocket that could create relevant neutron fluxes. However, the yields of some neutron-heavy 
isotopes have been shown in this paper to be very sensitive to the s-process in EMP stars. In 
these models the neutrons come almost exclusively from the He-shell flash convection zone 
and ^^Ne as the neutron source. In the future the role of the ^^C pocket for the s-process 
in general and for the EMP yields of species like the heavy magnesium isotopes must be 
studied in more detail. 

F. H. appreciates generous support from D. A. VandenBerg through his Operating Grant 
from the Natural Science and Engineering Research Council of Canada. This work would 
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13 


2 


104 


1 


C 


14 


( 00000, 00000) 


1 


N 


14 


c 


105 


1 


NI 


59 


( 00000, 00000) 


1 


CO 


59 


d 


C-burning 


59 


2 


C 


12 


( 00000, 1 PROT ) 


1 


NA 


23 


3 


60 


2 


C 


12 


( 00000, 1 HE 4) 


1 


NE 


20 


3 



1: NACRE adotped; 2: Horiguchi et al. (1996); 3: Caughlan & Fowler (1988); 4: Iliadis 
et al. (2001); 5: Bao et al. (2000a); 6: Beer et al. (2001); 7: Bao et al. (2000b); 8: 
Hauser-Feshbach (Jorissen & Goriely 2001); 9: see Jorissen & Goriely (2001); a: Wiescher 
et al. (1990); b: Koehler & O'Brien (1989); c: Takahashi & Yokoi (1987); d: Goriely 
(1999); e: Hauser-Feshbach (Jorissen & Goriely 2001) 
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Table 2. Model properties I 



ID 


Mini 


rric"' 


logL*'' 








AMDmax"^ 




tip' Miost'* 




Mo 


rric 


Lo 


Ro 


Ntd 


10^ yr 


10-2 Mo 


10-2 Mo 


yr 


Mo 


E82 


2 


0.626 ± 0.006 


3.96 


180 


6 


799.1 


1.171 


5.038 


82731 


1.34 


E84 


3 


0.814 ±0.006 


4.25 


250 


9 


277.1 


0.507 


3.465 


14116 


2.13 


E85 


4 


0.887 ±0.002 


4.40 


300 


10 


151.0 


0.367 


2.707 


8457 


3.05 


E79 


5 


0.947 ±0.001 


4.63 


360 


12 


97.2 


0.308 


2.422 


5281 


3.99 


E86 


6 


1.033 ±0.001 


4.70 


420 


15 


68.6 


0.140 


1.241 


2416 


4.83 


"mass 


of H-free core, range 


gives total ran^ 


le of core mass 


evolution, ^ 


'approximate average 





for the entire sequence, ^number of TPs with DUP, time at first TP, maximum 
dredged-up mass after a single TP, ^ total dredged-up mass all TPs, ^ average interpulse 
duration of TPs with DUP, ^ total mass lost=M* at first TP minus mnTP at end of AGB. 



Table 3. Model properties II 



TP txp" Tfbot'^ 2hes'^ ^hs'^ ^ceb"^ ttifbot^ ?^htp^ 'tt.dmax'* 



E82, M=2 M, 



1 


3.97£; - 02 


8.32 


8.10 


7.83 


6.31 


0.60550 


0.62559 


0.62740 


1.97251 


2 


8.23£; + 04 


8.42 


8.09 


7.81 


6.48 


0.60264 


0.62884 


0.62524 


1.96576 


3 


1.74£; + 05 


8.44 


7.93 


7.80 


6.65 


0.60278 


0.63041 


0.62403 


1.93878 


4 


2.70^ + 05 


8.46 


7.92 


7.80 


6.74 


0.60424 


0.63132 


0.62279 


1.86456 


5 


3.71^ + 05 


8.48 


7.93 


7.79 


6.79 


0.60230 


0.63193 


0.62088 


1.70939 


6 


4.76^ + 05 


8.49 


7.93 


7.79 


6.70 


0.60210 


0.63144 


0.61973 


1.44626 


7 


5.79^ + 05 


8.50 


7.94 


7.78 


6.20 


0.60065 


0.63111 


0.62021 


1.00773 


End of AGB 












0.62768 




0.62823 










E84, M=3 M, 












1 


2.05E - 04 


8.38 


8.16 


7.91 


6.70 


0.81151 


0.; 


31896 


0.81825 


2.94292 


2 


1.13^ + 04 


8.43 


8.14 


7.90 


6.99 


0.81156 


0.81945 


0.81711 


2.92614 


3 


2.40£; + 04 


8.46 


7.98 


7.89 


7.19 


0.81103 


0.81898 


0.81573 


2.88628 


4 


3.86£; + 04 


8.48 


7.99 


7.88 


7.34 


0.81034 


0.81819 


0.81430 


2.81285 


5 


5.46£; + 04 


8.50 


8.00 


7.88 


7.47 


0.80858 


0.81726 


0.81274 


2.68487 


6 


7.17^ + 04 


8.51 


8.01 


7.88 


7.56 


0.80752 


O.i 


31618 


0.81129 


2.50025 


7 


8.94E + 04 


8.52 


8.01 


7.88 


7.57 


0.80644 


O.i 


31508 


0.81011 


2.23381 


8 


1.08E + 05 


8.53 


8.08 


7.88 


7.41 


0.80549 


0.; 


31414 


0.80907 


1.87774 


9 


1.27E + 05 


8.54 


8.08 


7.87 


6.49 


0.80465 


o.i 


31331 


0.80830 


1.40821 


End of AGB 












o.i 


31251 




0.81343 










E85, M=4 M, 












1 


0.00^ + 00 


8.47 


8.35 


7.93 


7.72 


0.88627 


0.; 


39069 


0.88702 


3.93490 


2 


1.19E + 04 


8.48 


8.34 


7.94 


7.82 


0.88628 


0.; 


39012 


0.88756 


3.81435 


3 


2.01£; + 04 


8.50 


8.36 


7.94 


7.87 


0.88570 


0.88974 


0.88648 


3.65361 
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Table 3 — Continued 



TP 




Tfbot'' 


Thes'^ 




TcEB^ 


f 








4 


3.00E + 04 


8.50 


8.06 


7.94 


7.89 


0.88547 


0.88926 


0.88611 


3.36429 


5 


3.98E + 04 


8.52 


8.06 


7.94 


7.89 


0.88511 


0.88888 


0.88563 


3.00263 


6 


4.98E + 04 


8.53 


8.06 


7.94 


7.87 


0.88478 


0.88847 


0.88566 


2.55257 


7 


5.85E + 04 


8.54 


8.07 


7.93 


7.81 


0.88428 


0.88824 


0.88518 


2.12662 


8 


6.73£; + 04 


8.55 


8.08 


7.92 


7.58 


0.88380 


0.88785 


0.88481 


1.72477 


9 


7.62£; + 04 


8.55 


8.11 


7.91 


6.42 


0.88337 


0.88752 


0.88621 


1.36279 


10 


8.44£; + 04 


8.55 


8.00 


7.90 


5.23 


0.88469 


0.88843 


0.88747 


0.98297 


End of 


AGB 












0.88783 




0.88916 


E79, M=5Mq 


1 - 


-8.io£; + oi 


8.36 


8.20 


7.99 


7.38 


0.94796 


0.95020 


0.95024 


4.93359 


2 


2.64^ + 03 


8.40 


8.18 


7.99 


7.50 


0.94837 


0.95075 


0.95013 


4.91927 


3 


5.32£^ + 03 


8.44 


8.08 


7.97 


7.96 


0.94819 


0.95083 


0.94774 


4.90782 


4 


1.22E + 04 


8.48 


8.09 


7.98 


7.94 


0.94727 


0.94991 


0.94701 


4.79325 


5 


1.92^ + 04 


8.50 


8.10 


7.98 


7.89 


0.94655 


0.94907 


0.94645 


4.56126 


6 


2.64^ + 04 


8.51 


8.10 


7.98 


7.96 


0.94605 


0.94843 


0.94597 


4.13737 


7 


3.35£; + 04 


8.52 


8.10 


7.98 


7.95 


0.94558 


0.94795 


0.94560 


3.54257 


8 


4.02^; + 04 


8.52 


8.10 


7.98 


7.94 


0.94529 


0.94754 


0.94555 


2.90971 


9 


4.62£; + 04 


8.53 


8.10 


7.97 


7.92 


0.94510 


0.94740 


0.94550 


2.38893 


10 


5.14E + 04 


8.54 


8.11 


7.96 


7.86 


0.94494 


0.94730 


0.94547 


2.01317 


11 


5.61E + 04 


8.55 


8.11 


7.95 


7.71 


0.94468 


0.94723 


0.94531 


1.72970 


12 


6.10E + 04 


8.56 


8.12 


7.94 


7.35 


0.94440 


0.94706 


0.94524 


1.49238 


13 


6.58^ + 04 


8.56 


8.13 


7.94 


5.96 


0.94431 


0.94691 


0.94614 


1.25506 


14 


7.03E + 04 


8.56 








0.94486 


0.94756 




1.01774 


End of 


AGB 












0.94496 




0.94613 



E86, M=6 M, 



1 



-2.29E - 01 



8.37 8.22 8.03 7.92 1.03266 1.03362 1.03365 5.86578 
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Table 3 — Continued 



TP 




Tfbot'' 








f 

rriFBOT 


mHTP^ 




MJ 





i.UO£/ + Uo 


o.oy 


O.Zi 


a n/i 

O.U4 


/ .yo 


1 


.03297 


1, 


.03397 


1 


.03391 


K ac;Qn7 
o.oooU / 


q 
o 




O.'ii 


o.Zl 


a nc; 


7 Qa 
/ .yo 


1 


03323 


1 


.03426 


1 


.03414 


1^ a^aQi 
o.oooyi 




q 1 7 p _L nq 


o.4o 


8 1/1 
0. 14 


8 r\K 

o.UO 


8 nn 

o.UU 


1 


.03348 


1 


.03450 


1 


.03323 


80001 

O.oZZZi 





O.D4i!/ + Uo 


s /I 
O.40 


o.io 


8 nc; 
o.UO 


8 no 
o.Uz 


1 


.03327 


1 


.03424 


1 


.03306 


O.oy i ( i 


D 


1 nn ?? _L D/i 

i.UU-C/ -|- U4 


S AO 

o.4y 


O. i i 


8 nc; 

o.UO 


8 n'? 
o.UO 


1 


.03290 


1 


.03403 


1 


.03263 


1 fiQ7n 
0. ioy ( u 


7 


1.36^ + 04 


8.50 


8.12 


8.04 


8.03 


1 


.03255 


1 


.03358 


1 


.03246 


4.45788 


8 


1.70E + 04 


8.53 


8.12 


8.04 


8.02 


1, 


.03226 


1, 


.03340 


1, 


.03236 


3.74248 


9 


1.99E + 04 


8.53 


8.11 


8.03 


8.02 


1 


.03215 


1 


.03328 


1 


.03226 


3.15401 


10 


2.25£; + 04 


8.54 


8.12 


8.02 


8.00 


1 


.03201 


1 


.03319 


1 


.03229 


2.64631 


11 


2.50£; + 04 


8.54 


8.12 


8.02 


7.97 


1 


.03201 


1 


.03311 


1 


.03221 


2.33477 


12 


2.72£; + 04 


8.55 


8.12 


8.01 


7.94 


1 


.03192 


1 


.03312 


1 


.03221 


2.09245 


13 


2.92E + 04 


8.55 


8.12 


8.00 


7.90 


1 


.03189 


1 


.03309 


1 


.03222 


1.88476 


14 


3.12E + 04 


8.56 


8.13 


7.99 


7.79 


1 


.03186 


1 


.03311 


1 


.03227 


1.71168 


15 


3.32^ + 04 


8.57 


8.16 


7.98 


7.63 


1 


.03180 


1 


.03312 


1 


.03281 


1.56168 


16 


3.50E + 04 


8.57 


8.16 


7.98 


7.29 


1 


.03222 


1, 


.03354 


1, 


.03327 


1.45783 


17 


3.69£; + 04 


8.57 






6.20 


1 


.03272 


1 


.03402 


1 


.03377 


1.29629 


End of AGB 














1.034962 






1.065924 



" time since first TP, i = yr at about the time of the He-luminosity maximum during the 
first TP; Largest temperature at the bottom of the flash-convection zone; Temperature 



in the He-burning shell during the interpulse phase following the thermal pulse; 
Temperature in the H- burning shell during the interpulse phase following the TP; 
Temperature at the convective envelope bottom during the interpulse phase following the 
TP; ^ Mass coordinate of the flash-driven convection zone bottom; ^ Mass coordinate of 
the H-free core at the time of the TP; Mass coordinate of the convective envelope bottom 
at the end of the DUP phase; * Stellar mass at the TP. 
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Table 4. Average mass fractions in material returned to the ISM 



No 


species 




^ini 




E0082 




E0084 




E0085 




E0079 




E0086 








solar-scaled 




2Mn 




3Mn 

(J J.VJ.(V) 




4Mn 




5Mn 




\J J.VJ.(V) 


2 


PROT 


7, 


.6966E-01 


7, 


.260E-01 


7 


.336E-01 


6, 


.890E-01 


6, 


,600E-01 


6, 


.431E-01 


5 


HE 


4 


2, 


.3025E-01 


2, 


.569E-01 


2 


.588E-01 


3 


.070E-01 


3. 


,375E-01 


3 


.558E-01 


9 


C 


12 


1, 


.4750E-05 


1, 


.369E-02 


6 


.090E-03 


8, 


.980E-04 


4. 


,587E-04 


2, 


.159E-04 


10 


C 


13 


1 


.7750E-07 


5, 


.379E-07 


6 


.206E-06 


8, 


.565E-05 


4. 


,372E-05 


3 


.025E-05 


11 


N 


14 


4, 


.5551E-06 


2, 


.065E-05 


1 


.742E-05 


2, 


.331E-03 


1. 


,680E-03 


8, 


.002E-04 


12 


N 


15 


1 


.7950E-08 


6, 


.114E-09 


4 


.338E-09 


6, 


.465E-08 


4. 


,431E-08 


2, 


.023E-08 


13 





16 


4, 


.1401E-05 


2, 


.858E-03 


1 


.239E-03 


6, 


.llOE-04 


3, 


,008E-04 


1, 


.006E-04 


14 





17 


1, 


.6750E-08 


7, 


.067E-07 


5 


.453E-07 


9, 


.978E-07 


8. 


,188E-07 


4, 


.392E-07 


15 





18 


9, 


.3502E-08 


4, 


.984E-08 


4 


.930E-08 


1, 


.877E-09 


4, 


,783E-10 


6, 


.399E-11 


17 


NE 


20 


7, 


.6501E-06 


1, 


.548E-05 


1 


.837E-05 


1, 


.454E-05 


1. 


,168E-05 


9, 


.676E-06 


18 


NE 


21 


1, 


.9500E-08 


1, 


.002E-06 


9 


.041E-07 


2, 


.460E-07 


1. 


,020E-07 


3, 


.139E-08 


19 


NE 


22 


6, 


.1501E-07 


2, 


.733E-04 


5 


.500E-05 


1 


.272E-05 


4. 


,035E-06 


6, 


.782E-07 


21 


NA 


23 


1 


.6550E-07 


3 


.004E-06 


1 


.737E-06 


8, 


.356E-06 


3. 


,645E-06 


4, 


.964E-07 


22 


MG 


24 


2, 


.5250E-06 


3 


.351E-06 


2 


.914E-06 


3 


.016E-06 


2. 


,873E-06 


2, 


.884E-07 


23 


MG 


25 


3 


.3251E-07 


4, 


.961E-06 


5 


.390E-06 


3 


.291E-06 


2. 


,825E-06 


3 


.794E-06 


24 


MG 


26 


3 


.8101E-07 


6 


.937E-06 


1 


.155E-05 


7, 


.589E-06 


4. 


,351E-06 


2, 


.657E-06 


25 


AL26G 


1, 


.3000E-26 


4, 


.531E-10 


6 


.980E-10 


1, 


.109E-08 


5. 


,930E-08 


3 


.028E-07 


26 


AL 


27 


2, 


.8550E-07 


3 


.639E-07 


4 


.798E-07 


4, 


.679E-07 


4. 


,608E-07 


6, 


.706E-07 


27 


SI 


28 


3 


.2051E-06 


3 


.241E-06 


3 


.323E-06 


3 


.343E-06 


3. 


,293E-06 


3. 


,423E-06 


28 


SI 


29 


1, 


.6800E-07 


1, 


.842E-07 


1, 


.955E-07 


2, 


.035E-07 


1. 


,881E-07 


1. 


,780E-07 


29 


SI 


30 


1 


.1550E-07 


1 


.244E-07 


1 


.353E-07 


1, 


.442E-07 


1. 


,320E-07 


1. 


,240E-07 


47 


G 


63 


1 


.7000E-08 


1 


.499E-07 


8 


.597E-08 


5, 


.737E-08 


4. 


,277E-08 


3. 


,272E-08 


48 


L 


1 


1 


.OOOOE-50 


1 


.496E-08 


1 


.799E-08 


1, 


.207E-08 


6. 


,800E-09 


3. 


,551E-09 



a 
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Table 5. Yields EDITOR: this table only in the electronic version 



No 


species 




E0082 
2M0 




E0084 
3 Mo 




E0085 

4Mr^ 
^ J-VJ.0 




E0079 

5 Mr^ 




E0086 

W J.VJ.0 


2 


PROT 


-5 


.955E- 


■02 


-7, 


.794E-02 


-2 


.493E-01 


-4 


.335E-01 


-6 


.156E-01 


5 


HE 


4 


3 


.633E- 


■02 


6, 


.174E-02 


2 


.372E-01 


4 


.238E-01 


6, 


.104E-01 


9 


C 


12 


1, 


.867E- 


■02 


1, 


.314E-02 


2, 


.731E-03 


1, 


.754E-03 


9, 


.783E-04 


10 


c 


13 


4 


.921E- 


■07 


1 


.304E-05 


2 


.643E-04 


1 


.720E-04 


1 


.462E-04 


11 


N 


14 


2 


.197E- 


■05 


2, 


.784E-05 


7, 


.193E-03 


6 


.619E-03 


3 


.869E-03 


12 


N 


15 


-1 


.616E- 


■08 


-2, 


.945E-08 


1 


.444E-07 


1 


.042E-07 


1 


.lllE-08 


13 





16 


3 


.846E- 


■03 


2, 


.590E-03 


1 


.761E-03 


1 


.025E-03 


2, 


.879E-04 


14 





17 


9 


.421E- 


■07 


1, 


.143E-06 


3 


.034E-06 


3 


.169E-06 


2, 


.054E-06 


15 





18 


-5 


.962E- 


■08 


-9, 


.562E-08 


-2 


.833E-07 


-3 


.676E-07 


-4 


.544E-07 


17 


NE 


20 


1 


.069E- 


■05 


2, 


.318E-05 


2 


.131E-05 


1 


.592E-05 


9, 


.851E-06 


18 


NE 


21 


1 


.342E- 


■06 


1, 


.914E-06 


7, 


.005E-07 


3 


.259E-07 


5, 


.785E-08 


19 


NE 


22 


3 


.724E- 


■04 


1 


.177E-04 


3 


.744E-05 


1 


.351E-05 


3 


.071E-07 


21 


NA 


23 


3 


.875E- 


■06 


3 


.400E-06 


2 


.532E-05 


1 


.375E-05 


1 


.609E-06 


22 


MG 


24 


1 


.127E- 


■06 


8, 


.415E-07 


1 


.519E-06 


1 


.375E-06 


-1 


.088E-05 


23 


MG 


25 


6 


.320E- 


■06 


1 


.094E-05 


9 


.149E-06 


9 


.848E-06 


1 


.683E-05 


24 


MG 


26 


8 


.952E- 


■06 


2, 


.417E-05 


2 


.229E-05 


1 


.569E-05 


1, 


.107E-05 


25 


AL26G 


6 


.187E- 


■10 


1, 


.510E-09 


3 


.430E-08 


2 


.343E-07 


1, 


.473E-06 


26 


AL 


27 


1 


.070E- 


■07 


4, 


.202E-07 


5 


.640E-07 


6 


.925E-07 


1, 


.873E-06 


27 


SI 


28 


4 


.875E- 


■08 


2, 


.542E-07 


4 


.257E-07 


3 


.458E-07 


1, 


.058E-06 


28 


SI 


29 


2 


.213E- 


■08 


5, 


.959E-08 


1, 


.098E-07 


7 


.953E-08 


4, 


.866E-08 


29 


SI 


30 


1 


.211E- 


■08 


4, 


.288E-08 


8 


.859E-08 


6 


.501E-08 


4, 


.115E-08 


47 


G 


63 


1 


.815E- 


■07 


1 


.492E-07 


1 


.248E-07 


1 


.018E-07 


7, 


.646E-08 


48 


L 


1 


2 


.043E- 


■08 


3 


.891E-08 


3 


.732E-08 


2 


.687E-08 


1 


.727E-08 



a 



-29- 



Table 6. Structure and abundance evolution of computed tracks" 



model 








logL 


H 


^He 




5029 
5814 


2.00 
1.99 


-7.18-105 
3.23 • 10^ 


3.67 
3.66 


3.21 
3.48 


7.56 • 10"^ 
7.55 • 10"^ 


2.44 • 10-1 
2.454-10-1 


7.12-10-6 ... 
7.00 - 10-6 . . . 



« The complete tables ASPS.E82, ASPS.E84, ASPS.E85, ASPS.E79 and ASPS.E86 are 
available electronically. They contain for each sequence several structural and abundance 
quantities at 70 - 150 times during the TP-AGB evolution. Specifically these tables 



contain: model number, stellar mass in M0, age (t = at first TP), logTcfj, logL/L© and 
the surface abundance evolution of the 23 species given in Table 4. The models in each 
sequence ahvc been selected to reconstruct the abundance evolution. All details of the 
stellar parameter evolution (for example the exact stellar luminosity evolution after the 
thermal pulses) are only available from the complete evolution sequence, which can be 

obtained from the author. 
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Fig. 1. — Trends of stellar parameters of TP-AGB stellar models for two initial masses and 
a range of metallicities. The symbols represent values that roughly average the variations of 
the stellar parameters as a function of the TP cycle as well as the evolution towards the tip 
of the AGB. 
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Fig. 2. — Average core masses and stellar luminosities of the 5 Z=0.0001 sequences, with a 
linear m^. — logL-fit. As a comparison the core-mass luminosity relation by Blocker (1993) 
for non-HBB models and two points for Z=0 TP-AGB models during the last computed 
interpulse phase by Chieffi et al. (2001) are shown. 
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Fig. 3. — Core mass evolution. All panels show the same mass range of Amr = 0.025 Mg 




Fig. 4. — HRD for the pre-AGB evolution. The numeric labels indicate the initial stellar 
mass of the tracks. Note that higher-mass cases do not posses a RGB evolution phase. 
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Fig. 5. — Overabundances in the envelope just before the first TP, caused by the second 
dredge-up, and in the 2 M0 case through the first DUP. Xini is the solar-scaled initial abun- 
dance. Lines connect isotopic abundances for the respective stellar mass of the elements 
C, N and O, and for the Ne-Na and Mg-Al groups. EDITOR: THIS FIGURE DOUBLE 
COLUMN PLEASE 




Fig. 6. — Overabundances of CNO and Ne-Al isotopes in material returned to ISM. X[^i is 
solar-scaled and thus for each species log(Xav(/)/Xini) = [I/Fe], where the square brackets 
have their usual meaning of the logarithmic abundance ratio relativ to the solar ratio. Lines 
connect isotopic abundances for the respective stellar mass of the elements C, N and O, and 
for the Ne-Na and Mg-Al groups. EDITOR: THIS FIGURE DOUBLE COLUMN PLEASE 
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Fig. 7. — Abundance evolution of the most abundant CNO isotopes for all sequences during 
the TP-AGB. 
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Fig. 8. — Abundance evolution of the less abundant CNO isotopes. 
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Fig. 9. — Abundance evolution of the stable Ne and Na isotopes for all sequences during the 
TP-AGB evolution. 
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Fig. 10. — Abundance evolution of the stable Mg and Al isotopes for all sequences during 
the TP-AGB evolution. 
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Fig. 11. — Average mass fractions in ejected material according to this study compared to 
(Ventura et al. 2002). Apart from a small difference in metallicity the main difference is in 
dredge-up efficiency as a result of mixing assumptions (see text for details). 



